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Abstract 

The Lesche stability condition for the Shannon entropy [B. Lesche, J. Stat. Phys. 
27, 419 (1982)], represents a fundamental test, for its experimental robustness, 
for systems obeying the Maxwell-Boltzmann statistical mechanics. Of course, this 
stability condition must be satisfied by any entropic functional candidate to generate 
non-conventional statistical mechanics. In the present effort we show that the re- 
entropy, recently introduced in literature [G. Kaniadakis, Phys. Rev. E 66, 056125 
(2002)], satisfies the Lesche stability condition. 
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It is widely accepted that the Boltzmann-Gibbs statistical mechanics, while 
succeeding to describe systems in a stationary state characterized by an ergod- 
icity consistent with the thermal equilibrium, fails to depict of the statistical 
properties of anomalous systems as found in surface growth, anomalous dif- 
fusion in porous media, gravitation, Levy flights, fractals, turbulence physics, 
and economics [1,2]. 

Generalized statistical mechanics have been developed with the purpose of 
dealing with anomalous systems by means of the same mathematical tools 
used in the conventional statistical mechanics. This can be accomplished, in 
the case of trace-form entropies, with the introduction of an appropriate de- 
formed logarithm. Examples of trace-form entropies are the Tsallis entropy 
[3] and the Abe entropy [4]. In [5,6] some general properties of the deformed 
logarithms, leading to a generalization of the Boltzmann-Gibbs distributions, 
have been studied. 

Clearly, there are also entropies which are not of trace-form, e.g. the Renyi 
entropy [7], the Landsberg-Vedral entropy [8] and the escort entropy [9]. In 
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these cases the definition of these generalized entropies does not necessitate 
the preliminary introduction of a deformed logarithm. 

Recently, a non- conventional statistical mechanics has been introduced in 
[10,11] which preserves the epistemological and thermodynamics structure of 
the standard Boltzmann-Shannon-Gibbs theory. The proposed theory is based 
on the following trace-form entropy (k B = 1) 

N 

S K (p) = ~J2 Pi ln {K } (pJ > (!) 

i=l 



being p = {p t , i — 1, • • • , N} a discrete probability distribution normalized 

as 

N 

Eft = 1 > ( 2 ) 

1=1 



and the /t-logarithm is defined by 



The ^-logarithm is a monotonically increasing and concave function for k e 
(— 1, 1), being c?ln {K} (x)/(ia; > and gP ln {K} (x)/<ix 2 < 0, and satisfies the 
relation 

in w w=-v>Q • (4) 



In the limit of k — > Eq. (3) reduces to the standard logarithm and conse- 
quently Eq. (1) converges to the well-known Boltzmann-Shannon-Gibbs en- 
tropy. 

The inverse function of the ^-logarithm, namely the ^-exponential, is defined 

as 



exp {K} (x) = (v 1 + k 2 x 2 + k xj . (5) 



It is a positive, monotonically increasing and convex function which reduces 
to the standard exponential for k — > and satisfies the relation 

exp w (x) exp {K} (-x) = 1 . (6) 



2 



After maximizing the K-entropy given by Eq. (1) with the constraints (2) and 



N 



5><p, = <£>, (7) 



1=1 



being E. the energy of the i-th level with probability of occupation p., we 
obtain the probability distribution 



(8) 



where (3 — 1/T and —/3fx are the Lagrange multipliers associated to the 
constraints given by Eq.s (2) and (7), while the constants a and A depend 
on the parameter k and are given respectively by a = [(1 — k)/(1 + k)} 1 / 2k 
and A = y/l — k 2 . 

We recall now briefly the physical motivations which justify the introduction 
of ^-statistical mechanics. For a more detailed discussion we refer the reader 
to [10]. 

First we remark that the ^-exponential and the K-logarithm have the two 
following mathematical properties 

exp {R} (x) exp {R} (y) = exp {R} (x © y) , (9) 
ln w M =ln {lt} (x)© \n {K} (y) , (10) 

where the K-sum x © y is defined as 

x@y = x\Jl + K 2 y 2 + y^l + K 2 x 2 . (11) 



The above properties in the limit k — > reduce to the well-known relations 
exp(x) exp(y) = exp(x + y) and \n(xy) = ln(x) +ln(y) being y = x + y. Of 
course the K-sum plays a central role in the construction of the theory. In Ref. 
[10,11] it has been shown that the algebra constructed starting from the K-sum 
induces a /t-analysis. It results that the ^-derivative admits as eigenstate the 
^-exponential just as the ordinary exponential emerges as eigenstate of the 
ordinary derivative. At this point is evident that the physical motivation of 
the theory is strictly related with the physical mechanism originating the k- 
deformation. In Ref. [10] it has been shown that the K-sum emerges naturally 
within the Einstein special relativity. More precisely the K-sum is related to 
the relativistic sum of the velocities 

C V 1 + V 2 /-.^X 

1 + v 1 v 2 /c 2 
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through the relation 

p(uj ®p(v 2 ) =p(v 1 ® v 2 ) , 



(13) 



with re — 1/mc, where 

= = , (14) 



c 2 



is the relativistic momentum of a particle of rest mass m. Only in the clas- 
sic limit c — > oo the parameter re approaches zero and both the s-sum and 
the relativist sum of the velocities reduce to the ordinary sum. Thus the re- 
deformation is originated from the finite value of light speed and results to be 
a purely relativistic effect. 

Clearly when we consider a statistical system of relativistic particles the re- 
sum continues to play an important role. Due to the finite values of the light 
speed c 7^ oo, any signal and then any information propagates with a finite 
velocity in the system and therefore it results re ^ 0. 

The cosmic rays represents the most important example of relativistic statis- 
tical system which manifestly violates the Boltzmann statistics. In Ref. [10] 
it has been shown that the re-statistics can predict very well the experimental 
cosmic ray spectrum. This is an important test for the theory because the 
cosmic rays spectrum has a very large extension (13 decades in energy and 33 
decades in flux). 

Let us consider now statistical systems (physical, natural, economical, etc.) 
in which is involved a limiting quantity, like the light speed in the relativis- 
tic particle system. For these systems where the information propagates with 
finite velocity, it is reasonable suppose that the re-deformation can appear, 
so that the re-statistics results to be the most appropriate theory to describe 
these systems. 

An important question is if the re-entropy can be associated to a physical 
system in an experimentally detectable state. The problem can be stated in 
the following way. Let 0({x}) be a physical observable, depending on a set 
of parameters {x} = {x v i — 1, • • • , N} with sup[C({a;})] its maximal value 
achieved while varying {x}. 0({x}) is stable if, under a small change of the 
parameters {x} — > {x'}, with \{x'} — {x}\ < 5, the relative discrepancy 

_ Q({x>}) - Q({x}) 

sup[0({x})] ' 1 ) 



can be reduced to zero. 

This question has been considered, for the first time, by Lesche [12] in 1982 
both for the Boltzmann-Shannon-Gibbs entropy and for the Renyi entropy, 
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showing that the former is stable while the latter is unstable. Successively, 
Abe [13] showed that the Tsallis entropy is also consistent with the stability 
condition, while the Landsberg-Vedral entropy does not satisfy the stability 
condition. In [14] it has been shown that also the escort entropy does not satisfy 
the experimental robustness criterium. Finally, in [6] it has been considered 
the problem of Lesche stability condition from an arbitrary distribution func- 
tion. 

We remark that the Lesche inequality represents a condition for the stability 
of a system under particular sorts of perturbations. Clearly, the Renyi entropy, 
the Landsberg-Vedral entropy as well as the escort entropy, which fail to sat- 
isfy the Lesche stability condition, are worthy of investigation because present 
many other interesting properties. 

The purpose of the present contribution is to show that the K-entropy defined 
by Eq. (1) satisfies the Lesche stability condition and thus it can represent a 
well defined physical observable. 



Let us begin by introducing the following auxiliary function 



N 

K(P, s ) = S [Pi - a eX P {K }( _S ) 



i=l 



+ 



(16) 



being s a real positive parameter. The function [x]+ is defined in terms of 
the Heaviside unit step function 6{x) through [x] + = x9(x) so that it results 
[x] + = x for x > and [x] + = for x < 0. 

We remark that the function A K (p, s) has been considered already in the 
literature. Originally, in Ref. [12] A K (p, s) is defined using the exp(— s) at the 
place of a ex P {K} ( — s )- 111 Ref- [13] the function A K (p, s) is defined starting from 
the Tsallis exponential exp (s). Very recently, in Ref. [15], it was considered a 
very large family of functions having the structure of (16) where at the place 
of a exp {K} (— s) appears an arbitrary positive real function f(s). 
First we observe that from the relation \[x]+ — [y]+\ < \x — y\ it follows that 
Vs holds the inequality 

Kip, S )-K(q, s) < \\p-ql , (17) 



being 



N 



\\p-q\\i =J2\Pi ■ 
i=i 



(18) 



Second, following the procedure described in [12,13] we have that from the 
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definition of A K (p, s) follows 



N N 

E [ft ~ a ex p w (-*)]+ < Eft > 



(19) 



i=l 



which implies A K (p, s) < 1, being -fjr} a probability distribution normalized 
to one. 

On the other hand, by making use of the relation X)i[^J+ > we have 



N 



E [ft -« ex P{K} (-s) 
i=i 



> 



E(p,-«exp (K} H) 

i=l 

[l -aJV exp {R} (-s) 



(20) 



and posing s > ln {K} (iV), being a < 1, we can drop the notation [•]+ in the 
r.h.s of Eq. (20) obtaining 



A K (p, s) > 1-aN exp {K} (-s) 



(21) 



Combining Eq.s (19) and (21) we obtain 
1-aN exp {K} (s) < A K (p, s) < 1 . 



(22) 



Of course, this condition remains true for any other discrete probability dis- 
tribution {gj. Then, after changing the sign of Eq. (22), we obtain 



1< -A K (q, s) < -1 + a N exp {K} (s) . 



(23) 



Summing Eq.s (22) and (23) follows 

— olN exp {K} (s) < A K (p, s) - A K (q, s) < a N exp {K} (s) , (24) 



and than 

A K (p, s) - A K {q, s)\ < a N exp {K} (-s), 



(25) 



holding for s > ln {K} (A^). Eq.s (17) and (25) express two very important 
properties of the function A K (p, s) which will be used in the following. 
Now we show that the entropy S K (p) can be expressed in terms of the function 
A K (p, s). In fact, from the definition (1) and taking into account the identity 
a K (1 + k) — A, we have 
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N 



where L = 
we obtain 



i=i 



2k 



A f 



ft 



1 + k \a J 1 — k \a 



Pi 



A 



A? 



i=i 



1 fp.\ R 1 



Pi 

1 + k V«y 1-k \a 

TV P > 



= x ll j V l ds + J x — \n {K} (x)dx , 

i=l i=l x 

\n {K} (p./a). After introducing the substitution s = 



(26) 



ln {K} ( a; ) 



jv y n y 

i=1 i=1 L t 



ds 



A J ds - J [p { - a exp {K} (-s)] , 

o i=1 L ; 



(27) 



and finally we can write the entropy S K (p) in terms of the function A K (p, s) 
as follows 



S K (p) = \ J [l-A K (p,s)]ds-l. 



(28) 



-l/A 



Using Eq. (28) the difference of /t-entropy for two probability distributions, 
namely p — {p^} and q — {q^ can be written 



s k {p)-sm 



OO 

J [A K (p, s) - A K (q, s)} ds 

l/A 

30 

<A J \A K (p, s) - A K (q, s) 



(29) 



-l/A 



After splitting the integration interval in two parts, namely [—l/A, +oo) = 
[-1/A, e\{J[e, +oo), Eq. (29) becomes 
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S K (P)-SM\ <A / \A K ( P , s)-A H (q, s 

-1/A 

oo 

+A J\A K (p, s)-A K (q, s)\ds , 



(30) 



where the splitting point £ will be defined below. Now, using inequality (17) 
in the first integral of Eq. (30) and inequality (25) in the second integral of 
Eq. (30), we obtain: 



where the function Q K (£) is defined by 

l oo 

G K (£)= J \p-qlds + NaJexp {K} (-s)ds. 

-1/A i 

After performing the integrations, Eq. (32) can be written as 



(31) 



(32) 



N 



l + K 



1 + K 



+ 



1-K 



1 - K 



> + \\p-ql [£ + 



A 



(33) 



Let us choose for the parameter i the value £q which minimizes the function 
G K (£), therefore after posing 



dGM 



di 



we obtain 



=e 



, 



(34) 



In 



{K} 



aN 

lb -slli. 



(35) 



The condition s > ln {K} (iV), holding for Eq. (25), now requires £q > ln {K} (N) 
and then from Eq. (35) we obtain 



lb - Slli < a 



(36) 



By posing £ = £ in the expression of G K {i), the inequality (31) becomes 
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SM~ S M < lb- elli 



1 ( \\p-q\\i \ , 1 ( \\p-q\\i 

2 I N I 2 I N 



-v. [WW 



(37) 



and can be written in the compact form 



S K (p)-S H (q) < \\p-q\l 



In 



N 



{K} 



IIP - glli. 



+ 1 



(38) 



At this point we observe that f(x) = —x ln {K} (:r) is a positive and increasing 
function in the interval [0, a}. Furthermore we pose \\p — < 5 < a, consis- 
tently with Eq. (36). 

By introducing the maximum value of S K (p) corresponding to the uniform 
distribution p = {p. = 1/JV, i — 1, • • • , N} 



(39) 



from Eq. (38) we can evaluate the entropy relative discrepancy as 



SM-SM 



S 1 



Cl — K 



< 



5 l+K 2kS 
+ 



l_iV~ 2K N 2k — 1 N K — N~ K 



(40) 



Now, as customary we evaluate the inequality given by Eq. (40) in the ther- 
modynamic limit iV — > oo and obtain 



SM-SM 



< 5 



1- K 



(41) 



Finally, by introducing in Eq. (41) the positive quantity e = 5 1 ' K ', we obtain 
the Lesche inequality 



^max 



< e 



(42) 



Clearly in Eq. (42), e is a continuous function of 5, approaching for 5 — > 
when — 1 < k < 1 and therefore the /t-entropy S K satisfies the stability 
condition. 
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